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ABSTRACT 



Context. Weak gravitational lensing (WL) has been established as one of the most promising probes of cosmology. 
So far, most studies have exploited the shear effect of WL leading to coherent distortions of galaxy shapes. WL also 
introduces coherent magnifications. 

Aims. We want to detect this cosmic magnification effect (coherent magnification by the large-scale structure of the 
Universe) in large samples of high-redshift galaxies selected from the Deep part of the Canada- France-Hawaii- Telescope 
Legacy Survey (CFHTLS). 

Methods. Lyman-break galaxies (LBGs) selected by their colours to be at z=2.5-5, are used as a background sample 
and are cross-correlated to foreground lens galaxies, which are selected by accurate photometric redshifts (photo-z's). 
The signals of LBGs in difi'erent magnitude bins are compared to predictions from WL theory. An optimally weighted 
correlation function is estimated by taking into account the slope of external LBG luminosity functions. 
Results. For the first time, we detect cosmic magnification in a sample of normal galaxies. These background sources 
are also the ones with the highest redshifts so far used for WL measurements. The amplitude and angular dependence 
of the cross-correlation functions agree well with theoretical expectations and the lensing signal is detected with high 
significance. Avoiding low-redshift ranges in the foreground samples which might contaminate the LBG samples we 
can make a measurement that is virtually free of systematics. In particular, we detect an anti-correlation between faint 
LBGs and foreground galaxies which cannot be caused by redshift overlap. 

Conclusions. Cross-correlating LBGs (and in future also photo-z selected galaxies) as background sources to well un- 
derstood foreground samples based on accurate photo-z's will become a powerful cosmological probe in future large 
imaging surveys. 



^ 1. Introduction 



Many studies have shown so far that images of faint back- 
ground galaxies are coherently distorted by the gravita- 
tional lensing effect of the large-sca le structure in the 
line of sight. This cosm ic shear effect (|Munshi et al.|[2008t 
iHoekstra fc JainI l2008l for recent reviews see) has been 
identified to be one of the most promising approaches 
to study the properties of dark matter and dark energy 
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(CNRS) of France, and the University of Hawaii. This work is 
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France-Hawaii Telescope Legacy Survey, a collaborative project 
of NRC and CNRS. 



()Albrecht et aLll2006t Feacock et al.ll2006D . It relies on the 
assumption that the ellipticities of background galaxies are 
intrinsically randomly distributed. Cosmic shear introduces 
tiny, coherent distortions to this random distribution, which 
can b e measured (see e.g. iBeniamin et al.l 120071 : iFu et al.l 
I2008D . 

Besides this shear effect of gravitational lensing the im- 
ages of background galaxies are also subject to magnifica- 
tion by the large-scale structure of the Universe. However, 
no simple assumption about the intrinsic distribution of the 
fluxes of the background galaxies can be made. Rather, this 
distribution has to be measured from the data by averaging 
over large areas of the sky. 

The magnification effect of gravitational lensing to first 
order depends on the convergence only, which is the pro- 
jected mass along the line-of-sight. This effect is geometri- 
cal in nature enlarging the solid angle behind masses. This 
leads to two distinct effects, a dilution of the source density 
and a magnification of their fiuxes since lensing conserves 
the surface brightness. Now it depends on the intrinsic dis- 
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tribution of the fluxes whether the angular sky positions of 
the background galaxies of a given observed flux are pos- 
itively or negatively cross-correlated to the angular posi- 
tions of the foreground masses. If there are many more faint 
background galaxies than bright ones, i.e. steep magnitude 
number counts, the magnification wins over the dilution 
and a positive angular cross-correlation is expected. In the 
case of shallow number counts the dilution becomes dom- 
inant and a negative angular cross-correlation is expected. 
This general effect of gravitational lensing is called magni- 
fication bias and in the case of the large-scale structure of 
the Universe as the lens it is referred to as cosmic magnifi- 
cation. 

Measurements of the cosmic magnification signal can 
be interpreted in a simil ar way as signals o f galaxy-galaxy 
lensing (see e.g. [B rainer d et al.lll99 6: Hoek stra et al l 2004 
iMandelbaum et al.li200a iParker et al.ll2007t ). Thus, cosmic 
magnification can be directly employed to study the dark 
matter environment of galaxies constraining the galaxy 
bias, the total matter density, and the normalisation of the 
dark matter power spectrum. 

The biggest problem in using magnification as a cosmo- 
logical probe is to cleanly separate the fore- and the back- 
ground samples in redshift. If there is some redshift overlap 
between the samples - i.e. if the lens and source galaxies 
are physically close - they see the same dark matter field 
and will hence cluster with respect to each other. Then 
the angular cross-correlation signal is not a pure lensing 
signal anymore but it is contaminated by physical cross- 
correlations. 

Due to this complication the effect has only been con- 
vincingly measured with hig h significance in the Sloan 
Digital Sky Survey (SDSS) bv lScranton et all l|2005D using 
optically selected q uasars as background sources (see also 
iMenard et al.l [200 9. for a different estimator of the cos- 
mic magnification signal with the same data). For the long 
and controversial history of such measurements of quasar- 
galaxy cross-correlation we refer the reader to the references 
within that paper. 

Here we present the first measurement of the same effect 
on normal galaxies which have a much higher surface den- 
sity on the sky than quasars and can thus lead to much more 
accurate results. The Lyman-break technique allows for 
the selection of clean samples of high-redshift star-form ing 
galaxies from optical data. In iHildebrandt et all (|2009f ) we 
presented the largest survey of these galaxies to date. More 
than 80 000 LBG candidates with redshifts z = 2.5 — 5 were 
selected from the data of the Deep part of the CFHTLS 
and their clustering properties were measured. The same 
samples are used here to detect cosmic magnification by 
cross-correlating them to foreground galaxies. 

In Sect. Owe review the theoretical framework which is 
necessary to interpret the measurements. The data analysis 
is covered in Sect. [3] and the results are presented in Sect.lD 
Conclusions and an outlook to future applications of this 
method are given in Sect. [H Throughout the paper we as- 



sume Hq ~ 70 



km 



i Mpc 

and we use AB magnitudes. 



0.3, = 0.7, and cts = 0.8, 



2. Theoretical framework 

Let Nq{> f) be the unlensed cumulative number counts of 
background galaxies with fluxes larger than /. For simplic- 
ity we assume that all background galaxies are at the same 



redshift. Foreground structures will lead to a magnification 
/i at a particular position on the sky. The lensed cumula- 
tive number c ounts are related to the unlense d ones in the 
following way (jBartelmann fc Schneid"a1l200lh : 



(1) 



where the first fi corresponds to the dilution of the sample 
due to the enlargement of the solid angle behind the lens 
and the second /i corresponds to the brightening due to the 
enlargement of the sources. 

We assume that Nq^J) can be approximated by a power- 
law: 

7Vo(>/) = Ar", (2) 

with A being the amplitude and a being the slope of the 
power-law. The lensed cumulative number counts then be- 



N{> /)=/i"-iiVo(> /), 



(3) 



Thus, it depends on the slope a of No{> /) whether the sur- 
face density of background sources is increased or decreased 
near lenses where /i > 1. In the weak-lensing regime fi is 
close to unity so that we can write /i = 1 + S/i with S/j, <^ 1 
and a Taylor expansion yields 



1 + (a - l)Sfj.. 



(4) 



Using magnitudes instead of fluxes (i.e. substituting 
in = — 2.51og(/) + const.) it can easily be shown that a 
is related to the differential magnitude numbercounts: 



dlogn(m) 
Z.5 ; — a , 



(5) 



with n(rn) being the number counts of galaxies with mag- 
nitudes in the interval [m, m + dm]. By measuring the loga- 
rithmic slope of the magnitude number counts we can pre- 
dict over- or under-densities of background galaxies induced 
by lensing. 

Under the assumption of a linear biasing factor b for 
the foreground galaxies, the angular cross-correlation be- 
tween these foreground lenses and the background sources, 
Ws\{0), is related to the angular cross-correlation between 
magnification and matter density contrast, 'Wfj_s{9) by: 



Wsiie) = ia~l)bw^sid)- 



(6) 



For this result we employed Eqs.[3]&|31 the definition of the 
cross-correlation function, and the assumption of a linear 
relation between matter- and galaxy-density. 

We calculate Was(O) as described in 
iBartelmann fc SchneideTI (|200l[ ): 



3Hq f2ui 



kdk 



Ps{k,x')MfK{x')ke], (7) 



with X being the comoving distance, /k being 
the comoving angular diameter distance, Ws{x) — 

dx'Gsix') ^ being the weight function of the 

sources, Gi/s being the normalised distance distribution of 
the lenses/sources, o being the scale factor, Pg being the 
matter power spectrum, and Jg being the Oth-order Bessel 
function of the first kind. 
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In order to measure the signal of Eq.[n]from data it is of 
utmost importance to cleanly separate the sources from the 
lenses. Otherwise physical cross-correlations will swamp the 
tiny signal and make an interpretation in the framework of 
weak gravitational lensing impossible because these phys- 
ical cross-correlation are typically larger than the cosmic 
magnification signal by an order of magnitude. 

3. Data analysis 

3.1. The dataset 

The data used in this study are taken from the CFHTLS- 
Deep Survey, an imaging survey with MEGACAM@CFHT 
in the filters ugriz in four independent fields of 1 square 
degree each. In the framework of the CARS (CFHTLS- 
Archive-Research Survey) project we have collected all pub- 
licly available data until July 21, 2008. The data reduction 
is carried out wit h the THELI imaging reduct ion pipeline 
Erben et al.ll2005.) and is desc r ibed i n detail in iErben et al.l 



20091 ) and lHildebrandt erall (|2009l ). 



Multi-colour catalogues are created with SExtractor 
(|Bertin fc Arnout s 1996) in dual-image mode from images 
convolved to the same seeing. Photo-z's for all objects 
are estimated with a modified version of the BPZ code 
(|Beni'tej|2000[ ) incl uding a correction fo r galactic extinction 
(with the maps of ISchl eeel et al]ll998D . a re-calibration of 
the photometric zeropoints and the t emplate set w ith the 
help of spectroscopic redshifts (see ICapakI [2004( ). and a 
modified prior. For deta ils on the catalogue creation see 
iHildebrandt et al.l (|2009t) . 
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Fig. 1. Simulated redshift distributions of the three 
dropout samples (w-dropouts: solid^ ij-dropouts: dotted, 
r-dropouts: dashed] arbitrarily normalised but with cor- 
rect relative fractions). T he simulations are ba s ed on tem- 
plates from the library o f iBruzual fc ChariotI (|l993l ). See 
IHildebrandt et al] (j2009D for a detailed description of the 
simulations. 



3.2. The LBG catalogues 

In lHildebrandt et al.1 (|2009D we describe how we select large 
samples of LBGs from these data. Simulations are set up 
to identify regions in two-colour-space where high-redshift 
sources can be selected with high efficiency and low con- 
tamination. In this way we select ^ 34 000 w-dropouts at 
z ~ 3.2, ~ 36 000 5-dropouts at z ^ 3.8, and ~ 10 000 
r-dropouts at z ~ 4.7. The faintest LBGs that can be se- 
lected in that way from these data have measured total 
magnitudes of r = 27.6, i = 27.8, and z = 27.8 for the w-, 
g-, and r-dropouts, respectively. 

The simulated redshift distributions are displayed in 
Fig. [1] These suggest that the u-dropout sample is essen- 
tially free of any low-z contamination Q, whereas the g- 
dropout sample is contaminated by a small fraction (~ 4%) 
of low-z galaxies with redshifts < z < 0.5, and the r- 
dropout sample is contaminated very slightly (^^ 2.5%) by 
galaxies with redshifts 0.5 < z < 1.0. One of the big advan- 
tages of using LBGs as background sources is that we know 
the redshifts of the possible contaminants. In the remainder 
of the paper we try to avoid these redshift regions, which 
might be affected by some small amount of contamination, 
in our foreground lens samples in order not to mix the lens- 
ing signal with a signal from physical cross-correlations. 

For a given LBG background sample we approximate 
Gs of Eq. |7| by a Dirac delta-function at the mean redshift 
of the LBGs. This approximation is valid because the co- 
moving distance - the quantity on which the lensing signal 



^ There might be some contamination from galaxies at z ~ 1.5 
but the fraction should be very small with ~ 1% according to 
our simulations. 



depends - does not change appreciably over the range where 
the LBG redshift distribution is different from zero, i.e. the 
distributions in comoving distance are very narrow for the 
LBGs (in contrast to the redshift distributions). 

As described in Sect. |2l the amplitude of the cosmic 
magnification signal in the angular cross-correlation func- 
tion of low- and high-z galaxies depends on the slope of 
the number counts of the background sample. The num- 
ber counts of the three samples and of the combined uhg- 
dropout sample are shown in Fig.j^l For fainter magnitudes 
incompleteness sets in. While this incompleteness does not 
bias the measurement of the cross-correlation function (to 
first order) , it prevents a measurement of the slope of the 
number counts at the faint end, which is necessary to carry 
out the theoretical predictions. 

However, the number counts are closely related to the 
luminosity function (LF). For a complete sample of galaxies 
in a thin redshift slice the two curves are related to each 
other in magnitude by the distance modulus and in ampli- 
tude by the volume normalisation. Thus, using an external 
LBG-LF that has been properly corrected for incomplete- 
ness one can predict the slope of the number counts. The 
volume normalisation does not play a role here since we are 
only interested in the slope. 

The LBG-LF, ti>(M), was precise ly me asured in several 
studi es at z '--^ 3 (ISteidel et al. 11999 : Sawic ki fc Thompsonl 
'2OO I), z ~ 4 (ISteidel et airfT9 99: Ouchi et al.' '2004); 
Gi avalisco et all _|2004t ISawicki fc Thompson 200l 



Yo shida et al.l 120061: 'Bouwcn s et al l l2007l). and 



Pwata et aLll2003ll2007, : , Ouchi et al.ll2004l : lGiavalisco et al.l 
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Fig. 2. Number counts of the three dropout samples and of the combined u&g-dropout sample. 



12004 lYoshida et all |2006[) FI Assuming a Dirac delta- 
function for the redshift distribution of the LBGs, the 
slope of the number counts of a complete sample of LBGs 
equals the slope of the luminosity function. The parameter 
a introduced in Sect. [2] can thus be expressed as 

a(m) = 2.5dlogn(m)/dm = 2.5dlog$(Af)/dAf , (8) 

with m — M + DM + K and DM being the distance mod- 
ulus and K being the X-correction. 

For the theoretical p redictions we choose d i fferen t LFs 
bvlSteidel et all (I1999D . ISawicki fc Thompson! ()2006D . and 
iBouwens et al.l ( 2007 ) spanning a range of faint-end slope s 
and parametrised in the way described bv lSchechteil (|1976D : 



$(M)dM = 0.4$*ln(10) 



2q0.4(M.-M) 

xexp[-100-4(^'^--^'^)]dM. (9) 



Note that a(m) of Eq. [S] approaches the value of the 
Schechter LF parameter a^p for faint magnitudes. 

The Schechter function parameters of these external 
LFs are listed in Table [TJ We fit our number counts with a 
fourth order polynomial in magnitude only for comparison 
to the external LFs and not for prediction of the magnifi- 
cation signal. In Fig. [3] we show the adopted values of a — 1 
as a function of LBG magnitude for the different samples 
and the different LF measurements in the literature in com- 
parison to the fi tted p olynomial. We use the z ~ 4 values 
of ISteidel et"all (|l999f) and lSawicki fc Thompsol] (I2006D to 
estima te a z 5 LF and the z ^ 4 values of lBouwens et al.l 
()2007t) to estimate the z ~ 3 LF assuming no evolution. 

There is good agreement between the measurement and 
the literature LFs at the bright end (the r-dropouts suf- 
fer from small number statistics at z ^ 24.5). However, 
our measurement of the faint-end-slope suffers clearly from 
incompleteness. Thus, whenever a — 1 is needed in the re- 
mainder of the paper we present results for these three dif- 
ferent external LFs and do not use the number count slopes 
measured on our data. 



^ There is, however, some discrepancy between these litera- 
ture measurements. A proper measurement of the LBG-LF from 
the CFHTLS-Deep data themselves is under way (van der Burg 
et al. in preparation). Certainly it is better to calibrate the data 
internally and such a LF estimate would make the external cal- 
ibrators become redundant. 



Table 1. ISchechteil (119761 ) function parameters for the ex- 
ternal LFs 



Reference 

z 


~ 3 






Stcidcl ct al. (1999) 




-21.00 


-1.60 


Sawicki & Thomoson ( 


2006) 


-20.90 


-1.43 


Bouwens et al. (2007)'' 




-20.98 


-1.73 


z 


^ 4 






Steidel et al. (1999) 




-21.20 


-1.60^ 


Sawicki & Thomoson f 


2006) 


-21.00 


-1.26" 


Bouwens et al. (2007) 




-20.98 


-1.73 


z 


~ 5 






Steidel et al. (1999)" 




-21.20 


-1.60" 


Sawicki & Thompson ( 


2006)'' 


-21.00 


-1.26 


Bouwens et al. (2007) 




-20.64 


-1.66 



" This qlf is the faint-end-slope of the Schechter LF. The 
a(m) introduced in Eq. [2] approaches this value for faint mag- 
nitudes. 

'' extrapolated from 2 ~ 4 

" fixed 



3.3. Foreground samples 

The for eground samples ar e selec ted with the help of photo- 
z's. In iHildebrandt et ahl (120091 ) we showed that we can 
reach an accuracy of CTz/(i+z) = 0.033 rejecting only 1.6% 
of outliers in the magnitude range 17 < i < 24 if we filter 
for objects with a BPZ ODDS parameter of ODDS > 0.9. 
This filter basically rejects objects with a bimodal redshift- 
probability function. 

Three foreground redshift intervals have been selected: 
z = [0.1,1.0], z = [0.5,1.4] and z = [0.1, 0.5] U[l-0, 1.4] for 
cross-correlation with the u-, g-, and r-dropouts, respec- 
tively. For the analytical predictions, each redshift section 
is fitted with a three-parameters (riQ, zq, (Tz) Gaussian func- 
tion: 



n{z) = TiQ exp 



3.4. Clustering measurement 



(10) 



In the following, we call the catalogue of background LBGs 
Di containing A'di galaxies and the one of the foreground 
galaxies D2 with N^y^ galaxies. For the measurement of the 
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Fig. 3. Adopted values of a — 1 as a function of LBG magnit ude for the four background s amples. The do t ted, d ashed, 
and dash-dot t ed lin es correspond to the slopes of the LFs of ISawicki fc ThompsonI (|2006D . [Bouwens et al.l ()2007f ). and 
ISteidel et alj (jl999[ ). respectively, while the solid line corresponds to the slope of the polynomial fitted to the number 
counts of Fig. [H 



cross-correlation function we create one large random cata- 
logue, called R in the following, with the same areal geom- 
etry as the data catalogues and containing objects. We 
measure the angular cross-correlation function with a mod- 
ified v ersion of the estimator proposed by Landv fc Szalav, 
(fT99l : 

(11) 



D1D2 -D1R-D2R , ^ 
= RR + 



with D1D2 being the number of low-z-high-z galaxy pairs in 
the angular range [6, 9+S9] normalised by A'^Oi-^Da : DiR- be- 
ing the number of pairs between catalogue Dj and the ran- 
dom catalogue in that angular range normalised by TVd; A^r, 
and RR being the number of pairs in the random cata- 
logue in that angular range normalised by N^. By choos- 
ing Nb, ^ (by at least a factor of ten compared to 
the largest foreground samples) the shot noise introduced 
by the random catalogue can be suppressed. We use 10^ 
random points for each field that reduce to ~ 7 x 10^ af- 
ter masking. The masks used for the masking of the data 
catalogues are identical to the ones used for masking the 
random catalogues. Halos of bright stars are masked out as 
well as low-S/N regions (e.g. the borders of the stack that 
have lower S/N due to dithering) and regions affected by 
diffraction spikes or asteroid tracks. For a deta iled overview 
of the masking routines we refer the reader to lErben et al.l 
()2009f ) . This conservative masking approach results in a loss 
of ~ 30% of the area but ensures a highly uniform dataset 
wit h a homogene o us det ection and selection efficiency. 

iMenard et al.l (|2003D showed that the signal-to-noise 
of cosmic magnification measurements can be optimally 
boosted if an appropriate weight of a — 1 is put on each 
background galaxy. In this way the sources are weighed ac- 
cording to the expectations from the LF. Bright LBGs that 
are expected to be positively cross-correlated to the low-z 
lenses because of the steep exponential part of the LF get 
a positive weight. Faint LBGs from the shallow part of the 
LF that are expected to be anti-corrclated get a negative 
weight. And intermediately bright LBGs from parts of the 
LF where a — 1 w are down-weighed. We modify the 
estimator in the following way: 

w..^ DrD2-DrR-(w)D2R+(w)RR 

^'^ RR ' ^ ' 

with D^D2 and D^R being weighted pair counts, i.e. re- 
flecting the average of the weights of the LBGs in the se- 



lected pairs rather than the pure normalised number of 
pairs, and (w) being the average weight of the LBGs in 
the whole Di catalogue. 

Wc estimate the cross-correlation function separately 
for each of the four independent fields and calculate the 
mean u;(0). Furthermore, we draw ten jack-knife samples 
from the catalogue of each field and estimate the correlation 
function for all 40 of these. In order to take the correlation 
of the errors of data points at different angular scales prop- 
erly into account, the covariance matrix is then estimated 
in the following way from these jack-knife samples: 



N 



N 



(13) 



4. Results 

4.1. Cross-correlations in different magnitude bins 

First we cross-correlate LBGs in different magnitude bins 
to appropriate (i.e. non-overlapping) low-z samples to see 
if the signal agrees with the predictions. 

For the w-dropouts we choose the low-z range 0.1 < 
2phot < 1-0, for the g-dropouts we choose 0.5 < Zphot < 1-4, 
and for the r-dropouts we choose all galaxies with either 
0.1 < Zphot < 0.5 or 1.0 < Zphot < 1-4 (essentially a dou- 
ble peaked distribution). These choices are motivated by 
the simulated LBG redshift distributions shown in Fig. [TJ 
We exclude the low-redshift ranges that potentially con- 
taminate the LBG samples. Redshift beyond z = 1.4 are 
not considered for the lenses because between z = 1.4 and 
z — 2.5 we cannot expect our photo-z's to perform very well 
due to the lack of infrared filters. Furthermore, we restrict 
ourselves to magnitudes of i < 24 for the foreground sam- 
ple since without a deeper spectroscopic survey we cannot 
safely predict how the rate of catastrophic photo-z outliers 
develops for fainter galaxies. We apply an ODDS cut of 
ODDS > 0.8 as a compromise between accuracy and den- 
sity of the lens samples. 

In Fig. m the cross-correlation functions between the 
different source samples in different magnitude bins and 
the lens samples are shown. Errors are estimated from 
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jack-knife resampling. The magnitude bins were chosen in 
such a way that there are cases with predicted positive 
((a — 1) > 0) and negative {{a — 1) < 0) amphtudes as 
well as cases with an amplitude close to zero. For compar- 
ison also the predictions based on the three different LFs 
are plotted by using Eq. [6] in combination with the aver- 
age weight of the LBGs in the particular magnitude bin, 
(«-!)• 

4.2. Interpretation of the observed signal 

There is good qualitative agreement between the measured 
cross-correlation functions and the predictions from cosmic 
magnification. Bright LBGs show a strong positive cross- 
correlation to the foreground lenses, intermediately bright 
LBGs show a signal close to zero, and the faintest ones are 
anti-correlated. Especially the latter observation is a very 
strong argument for the lensing nature of the signal, since a 
physical cross-correlation caused by redshift overlap could 
not produce such an anti-correlation. 

The magnification predictions are pe rformed from Eq.[3 
using the non-linear power spectrum in iPeacock fc Dodda 
(|1996() with a biasing 6 = 1. There is a general tendency of 
an underestimation of the signal by the theoretical predic- 
tions. This can have different reasons: 

1. Lensing predictions are performed using the weak lens- 
ing approximation, i.e (k, 7) <§; 1 . Next order correc - 
tions are only of the order of 10% (jMenard et al.ll2003D . 

2. The power spectrum normalisation erg is probably also 
another minor contribution since it is unlikely that the 
real normalisation is very different from the fiducial 
0-8 = 0.8. 

3. The most likely explanation lies in the biasing param- 
eter b of the foreground galaxy population. Further in- 
vestigation will be necessary, in particular the com- 
bined analysis with the foreground auto-correlation 
function could remove any direct dependence on b 
(|Van Wae rbekc 2009). 

The predictions based on the LF measurements by 
ISawicki fc ThompsonI (I2006D seem to agr ee best with our 
data. The LFs bv iBouwens et al.l (|2007D estimated from 
space-based data with their very steep faint-end slopes 
do not yield the negative amplitudes observed in our 
cross-correlation functions of the faintest LBGs. However, 
we cut the LBG sam ples in observed magnitudes while 
IBouwens et al.l ()2007() account for the asymmetric scatter 
at the faint end introduced by rnagnitude errors called 
Eddington bias (|Eddingtonl Il913l : iTeerikorpil [2Q0l . The 
LBGs at the faint end of our samples have intrinsic magni- 
tudes that are on average fainter than the observed ones 
due to this effect. Taking this into account would lead 
to more negative values for (a — 1) and theoretical pre- 
dictions with a larger nega tive amplitude. Interestingly, 
ISawicki fc ThompsonI ()2006t ) do not correct for that asym- 
metric scatter. We suspect that the better agreement of 
our data with the predictions based on their LFs originates 
from that fact. 

The pred i ctions based on the LF estimates from 
Stcidel et al. (1999") lie in betw een the pr edictions from 
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It has also been reported by 
that cosmic variance can lead to a change in the shape of 
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Fig. 6. Correlation matrix (normalised covariance matrix) 
of the optimally-weighted cross-correlation function be- 
tween u-dropouts and galaxies with 0.1 < Zphot < 1-0. We 
display the scales 0'.3 < 9 < 15' used for the estimation of 
the total significance. 



the LF, especially at the faint end. Pencil-beam surveys 
in under dense fields tend to yield steeper slopes than the 
cosmic average. That may well be another explanation for 
the discrepancy between our results and the predictions 
based on the HST measurements. 

4.3. Optimally weighted cross-correlation functions 

Next, we estimate optim ally weighted c ross-c orrelation 
functioi is as introduced by [Menard et al.l (|2003l ) and also 
used in IScranton et al .' (2005). We weigh each galaxy with 
the a — 1 value corresponding to its magnitude and estimate 
the correlation function according to Eq. [T^l This is done 
three times fo r the three different sets of LFs. The results 
for the LFs bv lSawicki fc ThompsonI (|2006t ) are displayed in 
Fig. [5] together with the theoretical predictions. These are 
computed with Eq.|6]by taking the average squared weight, 
((a — 1)^), as the pre-factor. 

There is a similar tendency of the theoretical predic- 
tions being slightly lower than the observed signals, most 
serious for the u-dropouts. The same reasons as discussed 
in Sect. 14.21 apply here. 

Table |2] summarises the results for the normal as well 
as the weighted correlation functions. For the optimally 
weighted cross-correlation functions we report the total sig- 
nificance of the detection as computed with the help of the 
covariance matrix. See Fig. |B| for an example of such a ma- 
trix. 



4.4. Tests for systematics 

In order to test for possible systematics we select stars 
from our catalogues via a cut in magnitude and half-light- 
radius and cross-correlate these to our LBG samples. The 
amplitudes of the normal cross-correlation functions are 
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Fig. 4. Cross-correlation functions between the dropouts at different redshifts and with different magnitudes and the 
different foreground ga laxy sample s. The red, gree n, an d blue lines correspond to the predictions based on the LF slopes of 
ISawicki fc ThompsonI (12006 ) , Stei del et al. (1X9991 ) . and Bouwens et al. (|2007| ) . respectively. For some background samples 
the predictions by Sawicki fc ThompsonI "(|20n(iD and lSteidel et all (|1999h ~ are virtually identical so that the red and green 
curve lie on top of each other. 



mostly consistent with zero in all magnitude bins and for 
all LBG redshifts. The optimally-weighted cross-correlation 
functions are all consistent with zero as well. 

Furthermore, we checked the influence of the choice of 
the foreground sample. We included galaxies with photo-z 
estimates in regions where we would expect some contami- 
nation of the LBG samples. For example, including galax- 
ies with Zphot < 0.5 into the foreground sample that is 
cross-correlated to the g-dropouts leads to a boost in the 
amplitudes. In particular, the anti-correlations, which were 



observed before when excluding this \ow-z range, vanish. 
The signal turns positive for the faintest g-dropouts. This is 
in clear contradiction to the predicted lensing signal which 
should be negative because of the shallow slope of the LF at 
the faint end. Similarly, the negative signal for the faintest 
r-dropouts turns positive if galaxies with 0.5 < Zphot < 1-0 
are included in the foreground sample. These excess signals 
can be explained by redshift overlap leading to physical 
cross-correlations between the small number of contami- 
nants of the LBG samples and the foreground galaxies. 
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Fig. 5. Optimally weighted cross-correlation function between the complete dropout sample s and the different foregro und 
galaxy samples. The solid line correspond to the predictions based on the LF slopes of iSawicki &: ThompsonI (|2006f ). 



But then, the fact that we do see negative cross- 
correlations of the expected amplitude and angular depen- 
dence when these problematic foreground redshifts are ex- 
cluded from the low-z sample is a strong argument for the 
robustness of the analysis. While a small fraction of low-z 
contaminants is probably still present in our background 
LBG samples these do not change the amplitude or the 
shape of the signal but only add noise since they do not 
carry a lensing signal. 

5. Conclusions 

For the first time we detect cosmic magnification in sam- 
ples of normal galaxies. With the help of the Lyman- break 
technique we select background samples of high surface den- 
sity and large lensing efficiency (due to their high redshifts) 
from data of the CFHTLS-Deep survey. We cross-correlate 
these LBGs to low-z foreground galaxies which we select 
by accurate photo-z's. The expected signals are estimated 
by taking external LBG-LF estimates from the literature. 
There is good agreement between the observed signals and 



the theoretically predicted ones in amplitude as well as 
in angular dependence. Some deviations can be explained 
by Eddington bias, the linearisation of the magnification, 
and uncertainties in the cosmological parameters. The LBG 
samples used here represent the highest redshift population 
that has been used in weak gravitational lensing so far. 

Having proven that cosmic magnification with normal 
galaxies works in practice we plan to apply this technique 
to large imaging surveys in the future. In contrast to cos- 
mic magnification measurements with QSOs, using galaxies 
as sources has the advantage of much higher source den- 
sities. Also compared to cosmic shear and galaxy-galaxy- 
lensing this technique might prove competitive since the 
number of source galaxies with accurate magnitude mea- 
surements and photo-z's (the only requirements for magni- 
fication measurements) considerably exceeds the number of 
sources with accurate shape measurements. Although mag- 
nification measurements are less powerful than shear mea- 
surements for a given sample of galaxies, the larger densities 
that can be reached with future ground-based surveys will 
make precision measurements of cosmic magnification very 
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Table 2. Basic quantities of the samples used in the cross- 
correlation analysis. The subscripts rt. sti and bo refer to 
the luminosity funct i on est imat es of Sawicki fc Thornpsonl 
()2006D . ISteidel et al.l (|1999D . and lBouwens et al.l (|2007f l. re- 
spectively. 



sample 


No. of LBGs 


("-1>ST ("-l>St 

It-dropouts 


(«-1)bo 


(JST_ 


^ <^ 


l7Bo 


23.0 <r < 24.0 


714 


2.03 


1.97 


2.14 








24.0 < r < 25.0 


7263 


0.55 


0.62 


0.77 








25.0 <r < 25.8 


15524 


-0.02 


0.10 


0.24 








25.8 < r 


10679 


-0.29 


-0.14 


-0.01 








23.5 < r (opt. weights) 


34058 


0.17'' 


0.17'' 


0.24'' 


10.21 


10.47 


9.93 


(;-dropouts 


23.5 < i < 24.5 


990 


1.28 


1.28 


1.78 








24.5 < i < 25.5 


7557 


0.13 


0.32 


0.62 








25.5 < i < 26.5 


18948 


-0.36 


-0.08 


0.12 








26.5 < i 


8686 


-0.56 


-0.25 


-0.08 








23.5 < i (opt. weights) 


36181 


0.21'' 


O.IO'' 


0.20* 


6.95 


7.08 


5.46 


r-dropouts 


24.0 < 2 < 25.0 


520 


1.02 


1.07 


2.12 








25.0 < z< 26.0 


3482 


0.04 


0.25 


0.74 








26.0 < z < 26.8 


4768 


-0.37 


-0.09 


0.18 








26.8 < 2 


1685 


-0.53 


-0.23 


-0.05 








23.5 < 2 (opt. weights) 


10471 


0.20'' 


0.12'' 


0.51'' 


6.91 


5.94 


4.08 


ti&t;-dropouts 


23.5 < i < 24.5 


5256 


1.27 


1.27 


1.60 








24.5 < i < 25.5 


23212 


0.18 


0.33 


0.54 








25.5 < i < 26.3 


26619 


-0.26 


-0.06 


0.11 








26.3 < i 


14358 


-0.46 


-0.23 


-0.08 








23.5 < i (opt. weights) 


69445 


0.24'' 


0.20'' 


0.33'' 


8.98 


7.24 


4.96 



" detection significance 

^ average squared weight, ((a — 1)^) 



attractive. A precise calibration of the LFs of the source 
galaxies and a robust correction for observational bias like 
the Eddington bias is, however, mandatory to reach that 
goal. 

The cosmological constraints derived from cosmic mag- 
nification measurements will be complementary to other 
probes such as cosmic shear because their dependence on 
redshift is slightly different. This can potentially lead to 
breaking degeneracies in cosmological parameters. See also 
the theoretical companion paper by Van Waerbekc (2009) 
about this topic. Furthermore, cosmic magnification de- 
pends on completely different systematics. Measuring the 
same cosmological quantity with e.g. cosmic shear and cos- 
mic magnification simultaneously will be an extremely im- 
portant consistency check. Thus, cosmic magnification and 
cosmic shear can become a powerful combination in unrav- 
elling the mysteries of dark matter and dark energy. 
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